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In this paper we first apply the general analysis described in our first paper to a binary mixture 
of cyclohexane and n-hexane. We use the square gradient model for the continuous description 
of a non-equilibrium surface and obtain numerical profiles of various thermodynamic quantities in 
- - - various stationary state conditions. In the second part of this paper we focus on the verification 

of local equilibrium of the surface as described with excess quantities. We give a definition of the 
. temperature and chemical potential difference for the surface and verify that these quantities are 

independent of the choice of the dividing surface. We verify that the non-equilibrium surface can 
be described in terms of Gibbs excess densities which are in good approximation equal to their 
C_) , equilibrium values at the temperature and chemical potential difference of the surface. 
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I. INTRODUCTION 



In a previous article [l|, referred to as paper I, we have established the general approach for the square gradient 
' I description of the interface between two phases in non-equilibrium mixtures. We considered phenomena Uke temper- 

. ature, density and mass fraction gradients; heat and diffusion fluxes as well as evaporation or condensation fluxes 
^ through the interface. Some profiles were given, without going into details of the numerical procedures used to obtain 
I t , them. In this paper we will do this. 

d • In the general description of the interface one uses contributions to the Helmholtz free energy density proportional to 
^ the square of the density and mass fraction gradients. These contribution imply that it is not possible to use continuous 
I ■ local equilibrium thermodynamics in the interface, i.e. to calculate the local values of the various thermodynamic 
' ^ I parameters in terms of the local density, mass fractions and temperature only. Rowlinson and Widom (see [2, page 
g 43]) use the name point thermodynamics for this to distinguish it from other quasi- or local thermodynamic treatments. 
^ Given the non-autonomous nature of the square gradient model, it is sensible to question whether a description in terms 
I— I. of excess variables along the lines given by Gibbs [|, can be autonomous. Gibbs' treatment, though only given for 
equilibrium systems, suggested such an assumption. This would imply that the surface is a separate thermodynamic 
phase. Bakker fij and Guggenheim page 45] made this assumption, the validity of which was subsequently disputed 
by Defay and Prigogine [6[. We refer to Rowlinson and Widom 0, page 33] for a discussion of this point. In the 
theory of non-equihbrium thermodynamics of surfaces 0, BBIISj Gibbs' description in terms of excess variables has 
Cf^ ■ been used. It is then assumed that the Gibbs' description of the surface in terms of excess variables is autonomous, 
' or in other words that one can use this property, which we will call local equilibrium of the surface, to describe the 
surface. In earlier work (111] coauthored by one of us this property was verified for one-component systems. It is 
y—i the main objective of this paper to verify this property for binary mixtures. For details about the extension of the 
00 square gradient model to non-equilibrium systems we refer to paper I. For figures of typical density, mass fraction 
and temperature profiles, and a discussion thereof, we also refer to the first paper. In this paper we focus on the 
^ [ properties of the excess variables. 

We consider a flat interface between a binary liquid and it's vapor with the normal n — (1, 0, 0) pointing from the 
vapor to the liquid. We assume that all fluxes and gradients gradients are in the x— direction. Due to this all variables 
^ , depend only on the x— coordinate. We assume the fluid to be non-viscous, so that the viscous pressure tensor iTcfs = 0. 
■ - ' In Sec. [IT] we will give all the equations, which are required for the determination of the profiles. We shall only 
consider stationary states in this paper. We choose the system such that the gas is on the left hand side and the liquid 
is on the right hand side. The 2:-axes is directed from left to right and the gravitational acceleration g is directed 
towards the liquid. In Sec. [Illj we describe numerical procedure which have been used to solve these equations. We 
give the results in Sec. |IVj . We then proceed to the second part of the article, the verification of local equilibrium 
of the surface. In Sec. [V] we introduce excesses and surface variables and discuss in general the meaning of local 
equilibrium of the surface. In Sec. |VI| we give the results of the verification procedure. Finally, in Sec. |VII| we give 
concluding remarks. 
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II. COMPLETE SET OF EQUATIONS. 

In order to calculate the various profiles in the non-equilibrium mixture a complete set of equations is required. 
This set is given by the hydrodynamic balance (conservation) equations, the thermodynamic equations of state and 
the phenomenological force-flux relations. These equations were given (derived) and discussed in paper I, to which 
we refer. In this paper we will only list them for the stationary case. 

A. Conservation equations. 

In a stationary state the conservation equations take the following form. The law of mass conservation is 

d , 

— ipy) =0 
ax 

(III) 

(Ji+pCv) =0 

ax 

where p ^ pi + P2 and v = {pi vi + p2 ^2)1 P are the mass density and the baryccntric velocity. Furthermore ^ = pi/ p 
is the mass fraction of the first component and Ji = p^ (vi — v) is the diffusion flux of the first component relative 
to the barycentric frame of reference. Momentum conservation is given by 

(py"^ +P + jxx) ^ pg (II-2) 
ax ^ ' 

where we call 7^/3 the thermodynamic tension tensor, which will be defined below. Furthermore p and p -I- 72,3; are the 
pressures parallel and perpendicular to the interface, for the planar interface under consideration. For curved surfaces 
see paper I. Energy conservation is given by 

:^Je=0 (II.3) 

ax 

where Je = Jq + pe\ -\- p\ is the total energy flux, Jq is the heat flux, and e = u + (v^/2 ~ gx) is the total specific 
energy and u is the specific internal energy. 

B. Thermodynamic equations. 

The gradient model, discussed in the first paper Q gives the following expressions for the specific Helmholtz energy 
/, the specific internal energy m, the parallel pressure p, the chemical potential difference ^ = p,i — /i2, the chemical 
potential of the second component p = p,2 and the xx-element of the tension tensor j^x- 

f{x) =/o(r,p,e) + /c(p,c,p',C') 

u{x) = /o(r, p,0-T -^MT, P, + ^(p, ^, p', a 

p{^) = (/o(r, p, + np. t P\ ~P^[p ^Mp^ p'^ 

^{^) = ^ (/o(T, p, + /c(p, p', O) ~-pj^{p ^^P^ t p'^ 

p{x) = ^ (p (/o(r, p, + /c(p, P', m - i^i^) ^'P^{p ^^(^' p'' ^')) 

-fxx{x) =2p/c(p, p', n 

Here /o(T, p, £,) = fo{T, c, () is the specific Helmholtz energy of the homogeneous phase, which can for instance 
be derived from the equation of state. Furthermore JC{p, S., p', ^') is the gradient contribution, where the primes 



(11.4) 
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indicate a derivative with respect to x. Since the equation of state is usually given in molar quantities, it is convenient 
to use them here as well. Thus, c ~ p/M is the molar concentration, C = f M/Mi is the molar fraction, where 
M = Ml M2/(A/i + ^ (M2 - Ml)) = Af2 + C (^^1 - M2) is the molar mass of the mixture, and A'h and M2 are molar 
masses of each component. 

C. The homogeneous Helmholtz energy 

This energy is given by the following equation 

r,{T, c, C) = -RT In (^-^^L^ (1 - BiC) c)) - A{T, C) c (II.5) 
where the de Broglie wavelength A and the characteristic sum over internal degrees of freedom w are respectively. 

A(r, c) = hNA 



2n 



MRT 



w(T, C) 



C / .... i-C 



wi fMi\^^^\ ( W2 f M:^^'^ 



(II.6) 



Expressions for the characteristic sums over internal degrees of freedom for each component, wi and W2, are given 
in paper I. In this paper they are assumed to be independent of the temperature and the molar fractions, i.e. just 
constant numbers. 

The mixing rules for A and B are 

A{T, C) = an C + 2 012 C (1 - C) + "22 (1 - 0' .tt 7^ 

B{0 - &iC + 62(1-0 ^ ' 

with aij = y/ai aj , where ai as well as bi is a coefficient of a pure component i. We will assume in this paper that aij 
are independent of temperature. 

D. The gradient contribution 

This contribution is given by the following general expression for a binary mixture 

/c(p, p', C) ^ ^ {^ppiP, p'^ + 2 ^^pdp^ P' e + ^idp^ r') (11.8) 

The coefficients Kpp, Kpj and can be expressed in the gradient coefficients Kp^p^ and Kp2P2 for components 1 and 
2 in the following way (see paper I for details) 

l^ppiPj — {^Plpl ~ 2KpjP2 + KP2P2) C "I" 2 (Kpjp2 — Kp^pj) ^ + Kp2P2 

^p^iP: — i'^pipi ~ 2Kpip2 + '*P2P2) + {'^piP2 ~ '^pip-i) P (II-9) 
^^ii^Pi = {'^p\p\ ~ 2Kpip2 + i^pip-i) P 

where we use the mixing rule Kp^pj = \Ji^pipi i^P2P2 similar to the one for coefficients aij for the cross coefficient. We 
will assume Kp^p^ to be independent of the densities in this paper. 

With the above mixing rules the gradient contribution can be written in the form 

/C(p, ^, p', i') ^ (11.10) 
2p 

where k = Kp2p2 and q = p(l + £^ C), where e™ = = (V'^pipi ^ \/'^p2P2)I \/'^p2P2- Some of the quantities from 
Eq. (|II.4[) can be rewritten as 

p(x) =po - K (^9'^ + 9'7") 

p(x) =Aio-'«g" ("-^^^ 
V'(a;) =-00 -Eft Kg" 
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where po, fJ-o and ^po are values of the corresponding quantities in the homogeneous phase, which are found from 
Eq. (|II.4p by setting /C = 0. For a one-component fluid q equals the density. For the two-component mixture q plays 
a similar role as the density for the one-component fluid. We shall therefore refer to q as the order parameter. 

The size of the coefficient depends on the nature of the components of the mixture. In an organic mixture like 
cyclohexane and n-hexane, a mixture we will study in more detail in this paper, the components are very similar and 
as a consequence Sk is small. The order parameter is then in good approximation equal to the density. When the 
components are very different may be large and q may become in good approximation equal to the density of one 
of the components (for instance, the surface tensions of acetone and carbon disulfide are 22.7 and 31.3 gs"'^. Thus the 
values of corresponding Kp-p. are not very close and in their mixture the value of may be compared to 1). 

E. Phenomenological equations. 

In paper I we derived the general expression for the entropy production of a mixture in the interfacial region. For 
a non-viscous binary mixture which has only gradients and fluxes in the x— direction it takes the following form 



d I d ^ , , 



The resulting linear force-flux relations are: 



d_l_ 

dx T 

dx T 



(II.f3) 

= Rlq Jq — RuJl 



The resistivity coefficients Rqq, Rii and Rqi — Riq will in general depend on the densities, their gradients as well as 
on the temperature, so they vary through the interface. Expressions for the resistivity profiles in the interfacial region 
are not available. We model them, using the bulk values as the limiting value away from the surface and the order 
parameter profile as a modulatory curve. 

Rqqix) = + {R'qq ~ Rlq) qo{x) + aqq{R'qq + R^q) q,(x) 

Rqlix) =Rl, + {R'q,-Rl,)qo{x)+aqi{R'q,+Rl,)q,{x) (II. f 4) 

Rii{x) = + {Ril - qoix) + an{Rii + Rfi) qi{x) 



where 

<lo{x) ^ qi{x) = I '''/'^,|' — (11.15) 



^i^) -€q . , Wix)\ 



max 



are modulatory curves for resistivity profiles. Here gf„ and gf„ are the equilibrium coexistence values of the order 



eq ^^^^ Heq 

parameter of the gas and liquid respectively. Furthermore \q'eq{'^)\lmax the maximum value of the squared equilibrium 
order parameter gradient. For each resistivity profile R^ and R^ are the equilibrium coexistence resistivities of the 
gas and liquid phase respectively. Coefficients a^g, a^i, an control the size of the gradient term, which gives peaks 
in the resistivity profiles in the interfacial region. Such a peak is observed in molecular dynamic simulations of one 
-component fluids ^1^. 

Limiting coefficients R^ (where h is either g or t) are related to measurable transport coefficients in the bulk phases: 
thermal conductivity A'', diffusion coefficient and Soret coefficient s^. In the description of transport in the hulk 
phases it is convenient to use measurable heat fluxes 

Jq - Jl - Jl{h\ - hi) (11.16) 

where /i^ is a specific enthalpy of component i in phase h. Furthermore we used that j| = — J^. In the bulk phases 
the entropy production then takes the following form: 

, d 1 1 dibr 
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where we have suppressed the superscript b for now. The subscript T oi ip indicates that the gradient is calculated 
keeping the temperature constant. Using Gibbs-Duhem relation in a homogeneous phase at a constant pressure one 
can show, that 

dx {l~Odx ^ ■ ' 

After introducing measurable transport coefRcients, the force- flux relations derived from Eq. (|II.17p can be written in 
a form used in 13] : 



tJjJL/ '-^S U/U/ 

(11.19) 

Comparing Eq. ((TlT3)) and Eq. (|IlT9|1 in the bulk region we find the bulk values R of corresponding resistivity 
coefficients 



Rqq 



1 Dp 



Rqi ^Riq {hi ~h2) + DsTpH^-0T] (11.20) 

1 fDp,,^ ,,2 



Rii l^-^ (hi - + D ST p^il - {hi - h2) + XT 

where tp^ — {dtp/dS,), L = {\D p/ijj^) — {DsxpS, (1 — C))^ T. All the quantities in Eq. pi.20[) are taken in the specified 
bulk, either gas or liquid. 



III. SOLUTION PROCEDURE. 



The numerical procedure is similar to the one, described in [iJL however it has some differences. We will describe 
the special features below. We use the Matlab procedure bvp4c [1^ to solve the stationary boundary value problem. 
It requires a reasonable initial guess and boundary conditions. We use the equilibrium profile as the initial guess. We 
use a box of width 80 nm with the grid containing of 81 equidistantly spread points. 



A. Equilibrium profile. 

It is easier to describe equilibrium properties of the mixture using molar quantities. Everywhere in this subsection 
we will do this. The superscript v indicates a molar quantity. The total molar concentration and molar fraction of 
the first component are denoted by c and C respectively. 

Equilibrium coexistence is determined by the following system of equations 

Meg ~ A'o (-^e? ' '^eq ' Ceq ) ~ f^o{'^eq^ ^eq' Ceq) 

Veq = 4'^{Teq, Cf,, CI,) = ^UTeq, <, C'eq) (ULl) 
Peg = Po{Teq,4gXeq) = PO (T'eg , C^g , Cf g ) 

where V'o = (^/o /^C) 7 Mo = /o + c (9/o /dc) — V'o C and po = (? (9/g / 9c) are homogeneous chemical potentials and 
pressure, cf,, C^°' and Cg,, Cfg are coexistence density and mass fractions of gas and liquid respectively. 

Having 6 equations Eq. pil.ip and 8 unknowns cf,, Cfgj Cg,, Ceg and "i/jg,, /ig,, Peg, ?eg, an equilibrium state for 
two-phase two component mixture contains two free parameters. Particularly, the temperature and the molar fraction 
of the liquid phase are experimentally a reasonable choice. We have found, however, that it is more convenient to 
control Tgg and ^Aeg in the calculations, f/'eg changes monotonically with Cf, or Cfg, and it is therefore a good measure 
for the composition. Since — p!{ ~ the value of gives the difference of the chemical potentials of two 
components. 
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To obtain the equilibrium profiles Ceq{x) and Ceqi^) one needs to solve a system of two differential equations 

(111.2) 

Ve, =Vo^(c,C)-e>^(«'')" 

where q" = c(l + e^^C) = 9^2, = kM^ and ej;; = (1 + e^) (M1/M2) - 1 where Mi and M2 are the molar 
masses of the components. This system of equations is, in fact, singular, since coefficients of the higher derivatives 
are proportional. Thus, we can derive one algebraic equation instead of a differential one. 



ro (c, - < (c, c) = - < /^e, (ni.3) 

The bvp4c procedure takes only differential equations, so we have to transform Eq. piI.Sp to a differential one. This 
can be done easily by taking derivative of both sides. After some transformations, we have the following equation set 



rv 



(III.4) 



where subscripts C or c mean partial derivative of the corresponding quantity with respect to C or c. This is the system 
of 3 first order differential equations for 3 variables (, (f and (g"^)', which requires 3 boundary conditions. One of 
them is Eq. (jIII.3|) taken on one of the boundaries, which simply determines the integration constant for the second 
differential equation. The other two are {<fy(xg) = and [ifyixi) = which indicate the fact, that box boundaries 
are in the homogeneous region. 

Numerical procedure allows any values of the variables. However, not all the values are allowed physically. For 
instance, mol fraction C is bounded in the interval (0; 1) and molar concentration is bounded in the interval (0; 
where B is given in Eq. (|II.7|) . In order to avoid out-of-range problems, we use the function which safely maps unit 
interval to real axes and vice versa: 

u2r(u) = arcsin(2u - 1) and r2u(r) = 0.5(1 + sin(r)) (III. 5) 

Particularly, the actual variables, we provide to bvp4c procedure are 

Yi =u2r(C) 

Y2 =u2r(g70 (III.6) 
Fa -u2r'(g7gj;^„ (g'')VO 

where u2r'(u,M') — u' j \Ju — v? is simply the derivative of u2r and = B^^ niax(l, jej^l) is the limiting value for (f . 



B. Non-equilibrium profile. 

Non-equilibrium conditions are implemented when we change temperature or pressure from their equilibrium values. 
This results in mass and heat fluxes through the interface. The amount of matter will then change in the gas and 
liquid phase. We will put the system in such conditions, that the total contents of the box is constant and equal to 
the equilibrium contents. It means, that if some amount of liquid has been evaporated, the same amount of gas is 
condensed externally and put back into the liquid phase. 

We introduce the overall mass mix) = dy p{y) and the mass of the 1st component m^(x) = dy p{y) £,{y) , 
which obey the following equations by definition 

m'{x) = pix) 

m^(x) — p{x) £^{x) ^ ■ ' 
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We introduce the the overall mass flux J„j, the mass flux of the 1st component J^, the energy flux Jg and the 
"pressure" flux Jp-. 



1 /j,„y 1 / 1 ,2\ \ (iiii 



Jm 


= pv 








^ Jm 


Je 


= Jq^ 










Jp 


= Pl- - 






p 



2 \ p J p \ 2 



From Subsec. |II A| one can see, that all these fluxes are constant. 
From Eq. pi.lip we obtain 



,, 1 / 1 ,2 

q \ 2 



Kq" = — (V'o - tp) 

^ K. 



(IIL9) 



since p^{x) — p{x) + nq''^ . As in Eq. (jIII.2[) we have a singular set which leads to the algebraic equation 

^o-V'-— fpo-Pi + ^g'") =0 (III.IO) 



q V" ' 2 

Taking derivative of this equation with respect to coordinate we obtain the expression for the first derivative of the 
fraction (we do not give the exact expression, since it's very complicated) 

e'=e'(9,'7',e,T,V', Jm,Jc, Je,Jp) (III.ll) 

As a consequence we have 7 unknown variables, g, q', m, m^, ^, T, and 4 unknown fluxes Jm, J^, Je, Jp- This 
requires 7 first order differential equations and 11 boundary conditions (7 of them determine integration constants of 
differential equations and 4 of them determine constant fluxes). As differential equations we use Eq. pil.7p . Eq. pi.l3[) . 
one of Eq. pil.9[) and Eq. pil.lip . As boundary conditions we use the following: The first boundary condition is 
Eq. pil.lOp taken on one of the boundaries, which simply determine the integration constant for Eq. piLlip . The 4 
other conditions control the overall content of the box (particularly, it is the same as in equilibrium) 

m(xg) = 

mdxg) =0 (III 12) 

m{xe) = rrieq 

Here rrieq and m^. eg are the equilibrium values of the total overall mass and the total mass of the 1st component in 
the whole box. The 2 more conditions are 



q'ixg) =0 
q'ixi) =0 



(III.13) 



which indicate the fact, that box boundaries are in the homogeneous region. In contrast to the equilibrium case, the 
density in the non-equilibrium homogeneous region may vary with coordinate, so q' differs from zero on the boundaries, 
and, in fact, they do. The value of the q' is however small in the homogeneous region, comparing to the value in the 
surface region, so we may neglect it and use such approximation. This will lead to the wrong profile behavior only in 
the small vicinity on the boundary, which we will exclude from the further analysis. The 4 last boundary conditions 
reflect the conditions of the probable experiment. For instance, we may control the temperatures on the both side of 
the box, the pressure of the vapor side and the fraction on the liquid side. 



T{xi) =T' 

P±{Xg) =P 



a 



(III.14) 
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To solve these equations numerically we use the same technics as in the equilibrium case. All the variables should 
be properly scaled in order to make them to be the same order of magnitude. This balances the numerical residual 
and gives better solution result. We use the following variables 





= u2r(C) 


Y2 


= u2r(g/qoo) 








= raj {x*q*) 




= m^/{x*q*) 


Ye 


^ T*/T 


Yi 


= {yj/T){T*/r) 



where goo = 9So/^-^2 and scahng parameters x* = xe - Xg, T* = Teq, tp* ^ i^eq, q* Peq/fpeq- 



IV. RESULTS FOR THE TEMPERATURE AND CHEMICAL POTENTIAL PROFILES. 



In this section we show some profiles, obtained with the help of the above procedure. We choose a mixture of 
hexane and cyclohexane and give some of their properties relevant for our calculation in the table below. We note, 
that among them only the molar masses have been measured. There are number of problems to obtain the values of 
other material properties. We determined, for instance, the two van der Waals coefficients of the pure phases using 
their critical temperatures and pressures"'^. As a consequence the critical volumes per mole found in our description for 
the pure components differs substantially from the experimental value. For the mixture the van der Waals coefficients 
were then found using the mixing rules. We use the values of the molar mass and the van der Waals coefficients given 
in Table [J: 



TABLE I: The molar mass and the van der Waals coefficients 



component 


M, X 10"'' kg/mol 


a, J m'^/mol'^ 


b, xW" m''/mol 


1 


84.162 


2.195 


14.13 


2 


86.178 


2.495 


17.52 



Transfer coefficients of homogeneous fluids depend on temperature and densities, while these dependencies are not 
always available. We use typical constant values of these coefficients at the conditions, close to the above equilibrium 
conditions. The values of the heat conductivity A are well-tabulated and we take them from We take the typical 
value of the diffusion coefficient D for a liquid mixture from [3l and use the argument from [l^ p. 279] to obtain 
the typical value of the diffusion coefficient for a gas mixture. Another argument from P-279] is used to get the 
typical value of the Soret coefficient st- We use the data given in Table [II] 



TABLE 11: Transfer coefficients 







A, W/(m K) 


D, m'^/s 


ST, 


1/K 


phase \ component 


1 


2 








gas 


0.0140 


0.0157 


3.876 X lO-'' 


10" 


-4 


liquid 


0.1130 


0.1090 


3.876 X 10"^ 


10" 


-4 



^ In this we follow the example of the Handbook of Chemistry and Physics II6II rather than refs. [TH . Il4ll where Tc and Vc were used. 




together with the "mixing" rules for the heat conductivity 

(IV.l) 



The values of the gradient coefficients are not available at all. One can determine them comparing the actual value 
of the surface tension of a pure fluid with the one, calculated with a given Kp.p^ . For given conditions the value of the 
surface tension of the mixture is about 0.027 N/m. We therefore choose k'^ to be equal 12 x 10~^* J m^/mol^ and 
ej^ = 0.01. This gives values of the surface tension about 0.03 N/m. 

The equilibrium properties of the system are calculated^ at Teq = 330 K and tp'^^ = 700 J/mol. This gives 
Peg = 376095 Pa, /i^^ = -57098 J/mol, cf, = 153.23 mol/m^, cf^ = 4898.26 mol/m^, Cf, = 0.5519 and Cf, = 0.5934. 

The mixture is then perturbed from equilibrium for the following three cases: 1) setting equal to 0.98, 0.99, 
1.01, 1.02 of Teq and keeping T^, and equal to their equilibrium values, see Fig. [T]; 2) setting p^ equal to 0.98, 
0.99, 1.01, 1.02 of Peg and keeping T^, and C,^ equal to their equilibrium values, see Fig. [2]; 3) setting equal to 
0.98, 0.99, 1.01, 1.02 of (^q keeping , and p^ equal to their equilibrium values, see Fig. [3]. 



^ From now on we will use specific quantities per unit of mole in tlie description. We will omit the superscript u in the following sections. 
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V. LOCAL EQUILIBRIUM OF THE SURFACE 



In equilibrium it is possible to describe the surface in terms of Gibbs excess quantities |3| . One can treat a system 
of coexisting liquid and vapor as a three-phase system: liquid and vapor bulk phases and the surface phase. The 
surface has thermodynamic properties. The temperature and chemical potentials have the same equilibrium value as 
in the rest of the system. Furthermore the thermodynamic state of the surface is given by excess concentrations and 
thermodynamic potentials. Following Gibbs we have for the surface 



J eq 
9eq 



Ml, 



eq '-2,eq 
eq '^2,eq 



^eg Seq 
■ "teq ~^ ^eq ^eq 



^J■l,eq Cl,eq + M2,eg C^,e, + 7, 



(V.l) 



eq 



eq ^l,eq 



■ fJ-2. 



eq ^2,eq 



The superscript s indicates here the surface and equal to the excess of corresponding quantities. In these relations 
the temperature and the chemical potentials, which are the same everywhere, are independent of the choice of the 
dividing surface. The excesses depend on the choice of the dividing surface, in a such way that the above relations 
are true for any choice of the dividing surface. 

It is our aim in this paper to show that the surface in a non-equilibrium liquid-vapor system can also be described 
as a separate thermodynamic phase using the Gibbs excess quantities. We will call this property the local equilibrium 
of the surface. The property of local equilibrium for the surface implies that it is possible to define all thermodynamic 
properties of a surface such that they have their equilibrium coexistence values for any choice of the dividing surface 
given the temperature of the surface T'' and the chemical potential difference V'''- For this purpose we will define the 
excess densities and develop a method to obtain and f/'^ independent of the choice of the dividing surface below. 
The analysis will be done using the numerical solution of the system in stationary non-equilibrium states. 



A. Defining the excess densities. 



The definition of the excesses consists of 3 steps: determining the phase boundaries, defining the specified dividing 
surface and, in particular, defining the excesses. 

To determine the phase boundaries we will use the order parameter q. We introduce a small parameter j3 and define 
the /3-dependent boundary, x^'* between the vapor and the surface by 

where q^osip^, T) is the equation of state's value (no gradient contributions) of q for pressure p^, mass fraction ^ 
and temperature T. The /3-dependent boundary, Xg" , between the surface and the liquid is defined in the same way. 
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The numerical procedure calculates profiles only at specified grid points, which we provide to the procedure. That 
means that x^'* and x^'* can only be situated at points of the grid. We choose their position to be the last bulk point 
of the grid where the left hand side of Eq. (jV.2[) does not exceed the right hand side. In our calculations we will 
choose j3 = 10~^ and use a grid of 81 points. 

We shall also choose bulk boundaries near the box boundary where, because of the finite size of the box, the 
behavior of the profiles might be uncharacteristic. To avoid this effect, we do not consider the first 5 points of each 
phase close to these boundaries when we calculate the properties in these phases. The 6th point we call and the 
76th point x^ . 

The bulk gas therefore ranges from x^ to x^^^ and the bulk liquid ranges from x^'* to x^ . The surface therefore 
ranges from x^'^ to x^'^. In order to define excess quantities properly we always choose conditions such, that the bulk 

widths, x^0^ — x^ and x^ — a;^'* are larger then the surface width a;!'^ — a; a'*. 

In order to determine excess densities we need to extrapolate bulk profiles into the interfacial region. In equilibrium 
extrapolated bulk profiles are constants which are equal to the coexisting values of the corresponding quantities. 
Non-equilibrium bulk profiles are not constant. We fit the bulk profile with a polynomial of order n}, — 2 and use 
this polynomial to extrapolate non-equilibrium bulk profiles into the interfacial region. This is done with the help 
of Matlab functions polyf it and polyval. It is important to realize that the extrapolation of the bulk profiles 
introduces a certain error depending on the choice of /3 and rib m particular for non-equilibrium systems. 

The distances between commonly used dividing surfaces, such as for instance the equimolar surface and the surface 
of tension, are very small. Thus, if there occurs an error in the determination of a dividing surface using a course 
grid, that would lead to inaccurate results. We therefore divide each interval of the course grid between x^''* and x^''* 
in 10"* subintervals. This surface grid is used for all operations related to the surface. Within the interfacial region we 
interpolate all the profiles (which were obtained by extrapolation from the bulk to the surface region using the course 
grid) from the course grid to the surface grid using a polynomial of order rig = 3 with the help of Matlab functions 
polyf it and polyval. 

We can now define the excess (j) of any density (j){x) as a function of a dividing surface a;* as 

^{x')= ( ^ dx\^{x)-^^{x)<^{x' -x)-^\x)<^(x-x')\ (V.3) 

where 0^ and 0^ are the extrapolated gas and liquid profiles and 6(t) is the Heaviside function. The density (/> is per 
unit of volume and is per unit of surface area. In our calculations integration is performed using the trapezoidal 
method by Matlab function trapz. 

We can now define different dividing surfaces. The equimolar dividing surface x'^ is defined by the equation 'c(xf ) — 0. 
Analogously we define equimolar surfaces with respect to component 1 and 2: ci(a;'^^) — and C2(a:'^^) = 0, and the 
equidensity surface x'' . The surface of tension x^ is defined from the equation — x/p\\(xP^ = 0. All the 

densities are given as arrays on a coordinate grid, but not as continuous functions. Thus, in order to find the solution 
of an equation ip(x) = we calculate the values Lpi = Lp{xi) for each point Xi within the surface region and find the 
minimum of it's absolute value, mini(|0i|). Because of the discrete nature of the argument this value may not be 
equal to zero, but it will be the closest to zero among all other coordinate points. So we will call this point the root 
of the equation iy9(a;) = 0. We use the fine surface grid in this procedure. 

If follows from Eq. (|V3|) that 



d(t){xs 
dxs 

which we will use later. 



(b\x,)-^sixs) (V.4) 



B. Defining the temperature and chemical potential difference 

An equilibrium two-phase two-component mixture has two free parameters, for instance the temperature T and the 
chemical potential difference ip = — fi2- Local equilibrium of a surface implies, that in non-equilibrium it should be 
possible to define the temperature T'^ and the chemical potential difference of the surface. As found in ref. [ll| 
for the surface temperature in the one-component system, both and i/)^ should be independent of the choice of the 
dividing surface. 

The equilibrium temperature and chemical potential difference determine all other equilibrium properties of the 
surface. Thus, there is a bijection from Tgq and ipeq to any other set of independent excess variables Xi ^q and X2^eq, 



12 



so that one can use them equally well in order to characterize a surface. In non-equilibrium the temperature and 
chemical potential difference vary through the interfacial region, but as Xi ^e and X2,ne are excesses, they characterize 
the whole surface. If a non-equilibrium surface is in local equilibrium, there should exist the same bijection. This 
implies that given two independent non-equilibrium excesses Xi^ne and X2,ne one can determine the temperature T'^ 
and the chemical potential ip'^ of the whole surface. Thus one can calculate equilibrium tables of Xi^eqiTeq, V'eg) and 
X2,eq{Teq, tpeq) for different values of Teq and "tpeq and then determine temperature and chemical potential of a surface 

as = Teq{Xi^„e, X2,ne) and = 1peq{Xl.ne, -^2,ne)- 

As we want the temperature and chemical potential difference to be independent of the position of the dividing 
surface, we shall use excesses which are also independent of the position of a dividing surface in equilibrium for Xi 
and X2 ■ For two component mixture these independent variables are the surface tension 7 and the relative adsorption 
ri2. If the number of components is more then 2, additional relative adsorptions should be used. 

These quantities are well defined for equilibrium, but not for non-equilibrium. So we will define them first. In 
equilibrium the surface tension is defined as minus the excess of the parallel pressure "feq = —p\\- Alternatively 
one often uses the integral of p± — p\\{x) = ^xxix) across the interface: 7eg = J dx^xxix). Both definitions are 
equivalent in equilibrium since p±^ is constant through the interface and ^xxix) is identically zero in the bulk phases. 
In non-equilibrium ^xxix) may differ from zero in the bulk regions, however, and this makes the second definition 
inappropriate. We will therefore define the non-equilibrium surface tension using the standard definition 

7(a:^) = -p||(x^) (V.5) 

This quantity differs from 'jxx by the term equal to p±, which is usually small compared to p\\. 

The relative adsorption is defined as Ti2^eq — ci,eg — C2^eq {c{ — cf eq)/i'^2 eq ~ '^2 eg) in equilibrium where cf 
and are coexistence concentrations of the corresponding components. Since these quantities are not constant in 
non-equilibrium, we cannot use this definition directly. One can however see from Eq. (|V.4I) , that both, in equilibrium 
and non-equilibrium Ci'(x'*) = c^(a;^) — c\{x^\ where the prime indicates a spatial derivative. Since in equilibrium 
c\{x^) — cl(x^) = cf — cf we can use the following definition 



both in equilibrium and non-equilibrium. 

If the system is in local equilibrium we may write: 

7(x^) =7,,(^^v^) 
ri2(a:^) =di{x')-d2{xn 



cl,,{T^,r)-cieq{T',r) (V.7) 
ci.eq{T'.r)-4.eq{T'^r) 



Substituting the expressions for ^i{xs) and ri2(a:;s) from Eq. (|V.5[) and Eq. (|V.6[) into Eq. (|V.7[) we obtain the following 
relations 



c{{x^) - c\{x^) _ cl,g{T^,r)-cUq{T^r) 

ciix^) - 4ix^) 4,eg(^^V'^) - 4,,(T^ 



(V.8) 



This gives the bijection equations to determine and tp'^ from the actual non-equilibrium variables. As the left 
hand sides in Eq. (|V.8[) are in good approximation independent of the position of the dividing surface, T'^ and ip'^ are 
similarly independent on this position. 



C. Defining local equilibrium 



The other quantities required for the Gibbs description of the non-equilibrium surfaces we define in the following 
way. The surface chemical potentials are the equilibrium coexistence values determined via the procedure discussed 
in Subscc. ^JUl 

Ml =Ml,eg(T^?/'') /VQ^ 
=M2,eg(T^^^■) ^^-^^ 
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We define the surface extensive properties as^ 

(t)%x') = ^{x") (V.IO) 

The local equilibrium of a surface should be established for any choice of a dividing surface. The results of the 
calculations for any particular choice of a dividing surface may not be representative since they may be different differ 
for another choice of a dividing surface. Thus, the property of local equilibrium should be established for all dividing 
surfaces together. 




-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 
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FIG. 4: Equilibrium excesses at T = 330 K and ip = 700 J/mol. The vertical lines indicate the x'^^ , x~' , x'^, x'^^ dividing surfaces 
from left to right. 

Consider the profile of an excess <j>{x'') as a function of position of a dividing surface x'' . It follows from Eq. (IV. 4|) that 
the slope of the excess profile is equal to the difference between the extrapolated values of a profile of thermodynamic 
quantity (j). In equilibrium these values are constant and equal to the coexistence values. Thus, equilibrium excess 
profiles are straight lines, as one can see on Fig. [3]. Non-equilibrium profiles in the bulk phases are not constant. 
We construct the extrapolated profiles using nf,th order polynomials with ni, — 2. Resulting non-equilibrium excesses 
are therefore polynomials of the order n;, -I- 1 = 3, according to Eq. (jV.4[) . These profiles, for the most extreme case 
of non-equilibrium perturbation — 1.02 Tg^, are shown in Fig. [5]. Even though these profiles are polynomials of 
the 3rd order they are very close to straight lines. As one can see from Fig. [6] the variation in the slope is about 1% 
through the whole surface. It indicates that this non-equilibrium "state" is very close to an equilibrium one. 

We therefore develop the procedure to relate the non-equilibrium state to an equilibrium one by comparing thermo- 
dynamic quantities in equilibrium and in non-equilibrium. The comparison performed in one particular point of the 
surface may not be sufficient because it may suffer from artefacts peculiar to this particular surface. Moreover, any 
comparison performed in a particular point does not speak for the whole. We therefore compare the non-equilibrium 
surface with an equilibrium one for all dividing surfaces together. We will use the least square sum method for this. 

Consider a non-equilibrium thermodynamic excess and a quantity (x'* ; T, which is a combination of 

excesses and may depend on (T, ip) as parameters. We introduce the following measures of the difference of and r" 

Sts ix^;T,^)^\f {xt ) - (xt ; T, V') I (V. 1 1) 



^ Note, that for some quantities this definition diff'ers from the one, used in [Till . We will come back to this point later. 
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and 

i £ sur face 

(V.12) 

where N is the number of surface points. 

We say that and are the same in the surface if the value of fXts (r, is neghgible compared to the typical 
value of either |t'*(a;'*)| or |r*(a;*; T, If two quantities and r" are the same in the above sense, we say that the 
non-equilibrium state of the surface is characterized by surface temperature ^Y^s^ and chemical potential difference 

<S't^r4^{a;«}' V'la-s}) minS'ts_^s(T, V') (V.13) 

Here superscript s indicates that we speak about surface quantities only (as everywhere in this paper) and subscript 
{a:*} indicates that T^^s-^ and V'lj.Sj are the parameters for all dividing surfaces together (in contrast to the values 

T^{x^) and il)^{x'') determined from Eq. (IV. SP for each particular dividing surface x^). 

These definitions are easy to illustrate in equilibrium. For instance, for /i*^ gibbsi-'^'''^'^T'^) = /^i,e9(^i '0) cf + 
fJ'2,eq{T,ip)c^^^g + T slg it foUows from Eq. (jV.ip . that hlg{x^;Teq,il}eq) = Kq,g2bbsi^''^'^eq,ipeq)- Furthermore 

feqix'-.Teq.'Ipeq) ^ K qibbsi^" ' ^eq , Ipeq) ■ ThuS Sh^ , {x" ; Teq , eq) = and SgS , {x" ; Teq , eq) 0. 

^ ^' ^ eq eq, gibbs eq' eq, gibbs 

It is also true that Sh^ [Teq.ipeq) = TcmiShs j^s {T,ip) and cr^s j^s {T,tp) < agS j^s {T.tp). 

eq eq, gibbs eq eq, gibbs eq eq, gibbs eq eq, gibbs 

According to the above definitions i) h^^ ^-^f,^ and h^^ are the same quantities, but hl^ ^j^^^ and are not the same; 
ii) the equilibrium state is characterized by {Teq, ■(/'eg); as it should be. While this analysis is trivial in equilibrium, it 
is not trivial in non-equilibrium. 

Note, that while in equilibrium the conditions St,r{x'^] Teq, ipeq) = and St,r(Teq, ipeq) = niin St.t{T, ^) are equivalent, 
in general it does not follow in non-equilibrium from Eq. (|V.13|) that 

St,rix';T,^) = yx' (V.14) 
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FIG. 6: Slopes of non-equilibrium excesses for the case of perturbing = 1.02 Teg. The vertical lines indicate the surface 
boundaries and dividing surfaces. 

Thus Eq. (jV.14[) is not a good measure of the equaUty of the quantities and states in non-equilibrium. We may 
therefore speak about the equahty of thermodynamic quantities as well as about the state T" and ■0* of the surface 
in non-equilibrium only in the least square sense, as it is given in Eq. (|V.13[) . 

Within establishing the local equilibrium property of a non-equilibrium surface we would like to verify the following 
properties: i) the existence of the unique temperature and chemical potential difference of a non-equilibrium 
surface; ii) the validity of the Eq. (jV.ip in non-equilibrium at the surface's and V'*; hi) the possibility to determine 
all the properties of a non-equilibrium surface from equilibrium tables at the surface's and tp'^ . We do this in the 
following section. 

VI. VERIFICATION OF LOCAL EQUILIBRIUM. 

We calculate the the equilibrium properties (coexistence data, such as the pressure or bulk densities, as well as 
various excesses) of the system for the range of temperatures T = {325, 326, . . . , 340} K and the range of chemical 
potentials ip = {400, 450, 1000} J/mol. The value of a thermodynamic quantity at any point {T,iIj), which is 
between these is interpolated using the Matlab procedures interp2 and griddata. 

A. Surface temperature and chemical potential difference 

As was mentioned, in equilibrium both 7 and are independent of the location of the dividing surface x^. 
Given the above definitions, Eq. (jV.Sp and Eq. (|V.6p . we can calculate these quantities for non-equilibrium states. 
Calculations show, that even though 7 and are not exactly independent on a;** away from equilibrium, the relative 
deviation is so small (about 0.004% for 7 and 4% for in the worst case), that one can consider these quantities to 
be independent of the position of the dividing surface. Thus one may use them in order to find the temperature, T*, 
and the chemical potential difference, ip''^ , of the surface in non-equilibrium states, which will be independent of the 
position of the dividing surface. 
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Using Eq. (jV.12p and Eq. (jV.Sp together with Eq. (jV.4[) we construct the following expressions 

S'^(T,V) = ^ [p||(a;'') -p||,e,(r,V')]^ 



(VI.l) 



where the prime indicates the derivative with respect to x'^ 




Si(T^ tp) (where i is either 7 or ri2) should reach the minimum at T? and ip'} s, . We note however, that neither 

Si{T,ip) have a minimum at a single point (T*, V'")- There is a whole generatrix curve of minima Ci(T, V') = so 
the plot of Si{T,ip) is a valley. One can see it on Fig. [7]. Every point of the generatrix curve is the minimum 
point of Si{T,tp) along the direction "perpendicular" to this generatrix. If or c{{x^) and C2{x^) represent 

the corresponding profiles for some equilibrium state {Teq,ipeq), then Si(T,^) = and the generatrix is constant. 
Since these profiles are non-equilibrium profiles, the generatrix is not exactly constant but very close to it. Thus 
Ci{T,ip) = dSi{T,il))/dw, where w is a direction in T-^ plane which is perpendicular to generatrix. In fact, one 
should be careful speaking about directions, since no metric is defined in the T-ip plane. Thus we cannot introduce 
Vti/j so that Ci{T,tp) = \VTti;Si{T,^fj)\. In fact, w can be any direction which does not coincide or does not almost 
coincide with the direction of generatrix. In practice we find that we can use w = T, while using w = gives less 
accurate results. Thus we determine the minima curve from the equation 

^ (VI.2) 

Thus one needs two quantities S~. and S'r in order to determine Tf si and ip'f s-, uniquely. The surface temperature 
and chemical potential difference T^^sy and ''P^^sy are determined from the intersection of two minimum curves of 5^ 
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and Sr 

= 

} } (VI.3) 

= 

We calculate the temperatures and the chemical potential differences for different non-equilibrium conditions. They 
are outlined in Tables |IIII|V] . The first row of each table, corresponding to {x*}, gives and ip" calculated from 
Eq. (|VI.3|) . The following rows give, corresponding to different particular dividing surfaces, gives T** and tp^ calculated 
from Eq. (fyTS]) . 



TABLE III: Surface temperatures (K) and chemical potentials (J/mol) 
for the case of perturbing 







= 1.02 Te, 




= 0.98 Te, 


surface 




r 


rj-iS 


r 


{^n 


331.831 


770.53 


328.129 


650.92 


a;<= 


331.823 


769.51 


328.124 


650.29 




331.828 


770.22 


328.123 


650.21 


x"^ 


331.814 


767.97 


328.127 


650.43 




331.838 


771.86 


328.121 


650.1 



dT 



dT 



TABLE IV: Surface temperatures (K) and chemical potentials (J/mol) 
for the case of perturbing 





p9 


= 1.02 Peg 


p9 


= 0.98 pe. 


surface 




r 


rj-iS 


r 




330.796 


683.87 


329.059 


696.52 


x" 


330.8 


684.68 


329.063 


697.22 


x-^ 


330.799 


684.44 


329.063 


697.12 


x"^ 


330.804 


685.19 


329.065 


697.46 




330.795 


683.93 


329.061 


696.87 



TABLE V: Surface temperatures (K) and chemical potentials (J/mol) 
for the case of perturbing 





c 


= 1.02Ce, 




= 0.98C% 


surface 


rj-iS 


r 


ij-iS 


r 




329.577 


559.32 


330.24 


812.86 


x" 


329.598 


562.44 


330.242 


813.16 




329.584 


560.5 


330.253 


814.67 


x"^ 


329.63 


566.83 


330.219 


809.99 




329.554 


556.3 


330.278 


817.99 



Note, that T** and ip'' may be different from the continuous values in the interfacial region. 



B. The non-equilibrium Gibbs surface 



In this section we would like to verify that the surface quantities defined by Eq. (IV.10|) satisfy Eq. (|VI.4p with T'' 
and ^p" determined by Eq. (fOjl and Eq. (fVO)) . 

^<j^Ubs{T^r) (VI.4) 
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namely, with the definition Eq. (|V.9p . 



= ^lcl+ Ai2 + 



(VI.5) 



where the right hand side is 4'giiii,s(T'^ , tp^) of the corresponding quantity. Eq. (jVI.Sp is the non-equilibrium analogon 
of equilibrium Eq. (|V.1[) . 

In order to analyze the measure of validity of Eq. (|VI.4p we construct the quantity 



i £ surface 



(VI.6) 



for each thermodynamic potential h, u, /, g. fi^ ^.^^^ g ives the relative error of the determination of the surface 
quantity (p'' using the Gibbs excesses relations Eq. ()VI.5|) for for all dividing surfaces together, while €4,^^^^^ gives this 
error for particular dividing surface. We build £-ci>gMs (^j ^) fo^' T = T, 



for the whole surface. We build £0,^66^, (^i "0) both for T = T. 



^^sy, ip = V'^j.s}. determined from Eq. (|VI.3P only 
V'r si determined for the whole surface and for 
T = T'^^x'^), tp = ^p'^lx") determined from Eq. (|V.8|) for particular dividing surface. The values of the corresponding 
errors are listed in Tables |VIIIVIII| in Subsec. [A 1| . and are found to be small. 

As one can see, there is a variation in the value of the error for the different dividing surfaces. It is caused by two 
reasons. The first reason for this is slight variation in and tjj^ from Tables IIIHTV] for different dividing surfaces. 
The variation of each excess potential corresponds to the variation of and ■0" through these surfaces. Thus so do 
the relative errors. 

Another factor which influences the value of these errors is the actual value of an excess at a given dividing surface. 
If it is close to zero, then in the expression for e the small value is in denominator and it gives the huge value for the 
error. Particularly, g'^{x'^) ~ both in equilibrium and in non-equilibrium which makes the row corresponding to g 
at x'^ be uninformative and one should not take into account these data. 

We emphasize however that the overall error fi^ .^^^ represents the whole surface and thus do not suffer from the 
fact that some quantity is negligible at some dividing surface. There are such points for each potential 0, however 
their contribution to the whole error is negligible itself. So we can see, that if the particular dividing surface is far 
from zero point of 0, £4,^^^^^ gives the good measure of the error. While if the particular dividing surface is close to 
zero point of (j), e^^^^^^ fails to measure the error. One can see from Fig. [5] that the relative error e^fi^.^^^ indeed rises 
enormously at x'^. Particularly because of this fact the definition of the excess quantities in [ll| was different from 
Eq. dVTOl) . 



Another possible test is to compare the absolute error [(^"(a;*) 



i^gibbs y 



;T,ip'')\ with the deviation a^^.^^^iT^ip) 



defined in Eq. (|V.12[) . The calculations show that for the particular dividing surfaces the former quantity does not 
exceed the lateral, which indicates that all the absolute errors are actually within the trust region. 



C. Equilibrium tables 



In this subsection we will verify the possibility to determine all the properties of a non-equilibrium surface from 
equilibrium tables at the surface's T'* and t/j'^ . The surface chemical potentials and are already defined as their 
equilibrium values by Eq. (jV.9|) . So in this section we will verify the relation 

0^ = (^^^^) (VI.7) 

As in Subsec. [VI B| we compare the actual excess of a thermodynamic potential with the corresponding equilibrium 
value at given temperature and chemical potential of the surface. Before we do this a note has to be made. 

Under non-equilibrium conditions the profile of a quantity 0* is shifted with respect to the equilibrium one. One can 
see this in the example for (j) = h in Fig. [9]. The reason for this are fluxes caused by non-equilibrium perturbation. 
The whole surface is therefore shifted. One can clearly see that comparing the positions of the particular dividing 
surfaces on Fig. [5] and Fig. [S]. So the direct comparison of the profiles should be done not in the observer's 
frame of reference (OFO, which is used in all other calculations), but in the surface's frame of reference (SFO). 
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FIG. 8: The relative error e^gii,i,s(^^;^{x''}''^{x''}) perturbing = 1.02 Teq. The vertical lines indicate the 

surface boundaries and a;'' dividing surfaces. 
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FIG. 9: Non-equilibrium profile of h" for the case of perturbing = 1.02 Teq (solid line) compared to the profiles of hlq 
calculated from the equilibrium tables for different Teq and ipeq (dotted lines). The vertical lines indicate the surface boundaries 
and dividing surfaces for non-equilibrium case. 



The SFO is simply shifted with respect to the OFO, depending on the rate of non-equilibrium perturbations. Zero 
of the SFO is chosen at the reference surface, which can be either the equimolar surface, or any other physically 
sensible surface. If is the position of this surface in OFO and </'ofo('^OFo) profile of in OFO, then 

<^SFo(3^SFo) = <^OFo(3^0Fo) = '/'oFola^SFO + is the profile of (j)" in SFO. 

We can now determine to which equilibrium state the non-equilibrium one should correspond. Consider the following 
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definitions of E^j^^^ and e^^^^^ which have the same meaning as in Eq. (IVI.6P : 



H 2 



E 

i £ surface 



r{xn-€,i^';T,i^) 



(VI.8) 



(t)''{x'') 



for each thermodynamic potential h, u, /, g. and are the non-equiUbrium and equihbrium positions of the 
reference surface in OFO. The set {xf } is the non-equiUbrium surface grid and is used for both profiles. Since the 
width of an equilibrium surface may be not the same as the non-equilibrium one, the summation may exceed the 
formal boundaries of the equilibrium surface. This is not a problem however, since the equilibrium profile is the 
line with constant slope everywhere, as well as beyond the formal boundaries. We don't shift the surface grid in the 
definition of e^^^j^,^ (x"; T, because in that notation x^ means the particular dividing surface, while xf means the 
point of the surface grid. 

The values of the corresponding errors are listed in Tables |IXIIXI| in Subsec. |A 2j . and are found to be small. 

As in Subsec. |VI B| we see, that for the equimolar surface the relative error in g is huge. There is the same reason 
for this, namely that g^^Xc) ~ both in equilibrium and in non-equilibrium. This again makes the row corresponding 
to g at x*^ be uninformative and one should not take into account these data. 



VII. CONCLUSIONS 

The article continues the general analysis started in ^ . Here we focus on a specified mixture and develop further 
its properties. We choose a binary mixture of cyclohexane and n-hexane as the system and describe in details how 
to implement the general analysis presented in We build numerical procedure for solving the resulted system of 
differential equations. The resulted profiles of continuous variables are presented in Sec. |IVj and in the first article. 
We see, in particular, that a two component mixture develop the teraperature profile in the surface region which is 
similar to the temperature profile obtained for one-component system [14| . Another characteristic of a binary mixture 
is the difference between chemical potentials of components. The behavior of the profile of ip in non-equilibrium 
steady-states is different. It has different values in the different bulk phases and we observe a transition from the one 
value to the other in the surface region. 

We then proceed to verify the local equilibrium of the surface. This property means that a surface under non- 
equilibrium steady-state conditions can be described as equilibrium one in terms of the Gibbs excess densities. We 
have discussed the meaning of the surface quantities in non-equilibrium and established the systematic procedure to 
obtain them. We were in particular focused on i) the existence of the unique surface's temperature and chemical 
potential difference; ii) the validity of the relations between thermodynamic Gibbs excesses in a non-equilibrium 
surface; iii) the correspondence between the non-equilibrium and equilibrium properties of the surface. It was possible 
to verify that one can speak about these statements independently on the choice of the dividing surface. Similar 
results for the one-component system were obtained in [ll[ . 

The extrapolations procedure is numerical, and contains therefore a certain error. We may not expect this error to 
be negligible, not only because of numerical inaccuracy, but also because of the non-equilibrium nature of the system. 
If the error is within a reasonable range, we will consider this as a satisfactory verification of local equilibrium. 

The main part of the analysis in the interfacial region is introducing the excesses of thermodynamic densities, which 
are constructed with the help of extrapolated bulk profiles. In contrast to equilibrium, non-equilibrium bulk profiles 
are not constants, and therefore their extrapolation to the surface region is not always accurate. The accuracy of 
extrapolation lowers when the surface width increases. Apparent small deviations from local equilibrium are therefore 
to some extent an artefact of the inaccuracy of the extrapolation. 

In the description of the surface excess densities it may happen that for particular choice of the dividing surface 
not one but several of the excesses are negligible. This increases the relative error enormously while the absolute 
error remains finite and more or less constant. In order to avoid this problem we consider the excesses for all dividing 
surfaces together, rather then for a particular dividing surface. Particularly, in the definition of the excess Gibbs 
energy was chosen differently because this excess was very small for the equimolar surface. We have shown in this 
paper why this is not needed. 

One can see from these data, that within different ways of perturbing mixture from equilibrium the biggest error 
comes when one perturbs the temperature on the liquid side. This is the most extreme condition for the mixture 
being in non-equilibrium. While the relative temperature perturbation is only 2%, the resulting temperature gradient 
is about 10* K/m which is very far beyond ordinary non-equilibrium conditions. The other perturbations make the 
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validity of local equilibrium for the surface more precise. Similarly smaller perturbations make the validity of local 
equilibrium also more precise. 

We therefore conclude that the local equilibrium of the surface is valid with a reasonable accuracy. 

APPENDIX A: EXCESSES' ERRORS 
1. Gibbs excesses' relative errors 



TABLE VI: Gibbs excesses relative error for the case of perturbing 
in percent 









T" = 1.02 Teg 






= 0.98 Teg 




error 




for T\x''),'4i\x') 






for T\x'),'i>\x'') 
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0.01328 


- 


0.033276 




- 






0.023799 


0.023419 


0.064597 




0.065002 






U.UzUdUo 


0.020737 


0.05d5o5 




A ri Tl AO 

0.057148 






0.035878 


0.033571 


0.085212 




0.085204 






0.015606 


0.016541 


0.047847 




0.048578 
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£ 


0.0071257 


- 


0.026714 




- 






0.016729 


0.016462 


0.045109 




0.045392 






0.015059 


0.015156 


0.041048 




0.041457 






0.022033 


0.020616 
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0.012161 
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0.20039 


- 
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e^{x ) 
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7.3666 


8.7959 
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£ 


0.36323 




0.15487 
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73.468 
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e^ix'^) 
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2.7242 
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1.4054 


1.3983 


1.9592 
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TABLE VII: ' 


Gibbs excesses relative error 


for the case of perturbing -p^ 






in percent 
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0.00083823 






0.0025175 
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0.0015513 
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TABLE VIII: Gibbs excesses relative 


error for the case of perturbing 






in percent 
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0.0017677 


0.0025842 


0.0027086 


f 


£ 


0.0063873 




0.0026805 








0.17101 
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2. 


Equilibrium table excesses' relative errors 








TABLE IX: Equilibrium table excesses relative error for the case 


of per- 






turbing T in percent 
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5.2046 

3.0361 


0.35744 
0.053575 
4.1104 
4.692 


0.3558 
0.051756 
4.1143 
4.6875 


f 


£ 
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TABLE X: 
turbing 


Equilibrium table excesses 
in percent 


relative error for the case 


of per- 
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0.047434 
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1.0589 
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TABLE XI; Equilibrium table excesses relative error for the case of per- 
turbing in percent 
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2.2203 
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